Field-induced phase transitions in a Kondo insulator 
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We study the magnetic-field effect on a Kondo insulator by exploiting the periodic Anderson model 
with the Zeeman term. The analysis using dynamical mean field theory combined with quantum 
Monte Carlo simulations determines the detailed phase diagram at finite temperatures. At low 
temperatures, the magnetic field drives the Kondo insulator to a transverse antiferromagnetic phase, 
which further enters a polarized metallic phase at higher fields. The antiferromagnetic transition 
temperature T c takes a maximum when the Zeeman energy is nearly equal to the quasi-particle 
gap. In the paramagnetic phase above T c , we find that the electron mass gets largest around the 
field where the quasi-particle gap is closed. It is also shown that the induced moment of conduction 
electrons changes its direction from antiparallel to parallel to the field. 

PACS numbers: 71.27.+a 71.10.Fd 71.30.+h 75.30.Mb 
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I. INTRODUCTION 

There has been a continued interest in a class of 
compounds called Kondo insulators or heavy fermion 
semiconductors, 1 which develop an insulating gap at low 
temperatures via strong correlation effects. CeRhAs, 
Ce3Bi4Pt3, YbBi2, and SmBg are well known examples, 
which possess small gaps of the order of 1-100 meV in spin 
and charge excitations. The gap formation in Kondo in- 
sulators is attributed to the renormalized hybridization 
between a broad conduction band and a nearly flat /- 
electron band with strong correlations. In some Kondo 
insulators, enhanced antiferromagnetic (AF) correlations 
among localized /-moments dominate the Kondo singlet 
formation, resulting in an AF long-range order & Systems 
located around such magnetic instability provide hot top- 
ics in quantum critical phenomenal 

The application of a magnetic field to the Kondo insu- 
lators also realizes quantum critical phenomena, which 
have been investigated extensively. Experiments on 
Ce3Bi4Pt3j' 1 YbBi2j§ and SmB^ 6 in high magnetic fields 
indicate closure of the Kondo-insulating gap, exemplify- 
ing a transition from the Kondo insulator to a correlated 
metal. If the Kondo insulator is in the proximity of mag- 
netic instability, the local singlet formation gets weak, 
and the applied field may possibly trigger a phase tran- 
sition to the AF ordered state before it becomes a metal. 

The periodic Anderson model (PAM) at half filling 
may be a simplified model to describe the Kondo insu- 
lating phased and the AF phase, depending on the c-f 
hybridization strength V and the /-/ Coulomb interac- 
tion U. The magnetic instability of the PAM has been 
investigated by means of various methods such as slave- 
boson mean field theory&2*i2i and dynamical mean field 
theory (DMFT). 11 ! 12 ! 13 ! 14 The effect of the magnetic field 
on the PAM has also been investigated) 15 : 16 : 17 ! 18 ! 19 ! 20 and 
the phase transition from the Kondo insulator to a para- 
magnetic metal, which is naively expected, has been dis- 
cussed. 

Recently, Beach et. alM> have studied the magnetic- 
field effect on the two-dimensional Kondo lattice model 



using large-iV mean field theory and quantum Monte 
Carlo (QMC) simulations. Also, Milat et. al. have pre- 
sented a mean field analysis of the PAM in the small 
U region and QMC simulations of the Kondo lattice 
model, 22 They have found that the magnetic field in- 
duces a second-order phase transition from the param- 
agnetic to the transverse AF phase in the Kondo insu- 
lator. Since these works have been concerned with the 
two-dimensional systems, it is desirable to extend the in- 
vestigation to three-dimensional cases. Also, a detailed 
study including the paramagnetic phase at finite temper- 
atures is expected to provide further interesting proper- 
ties in the Kondo insulator in a magnetic field. 

In this paper, we study field-induced AF phase tran- 
sitions of the Kondo insulator in the whole temperature 
regime by using the PAM at half filling. By exploiting 
DMFT combined with the QMC method^ we show that 
a magnetic field induces a transverse AF order in the 
Kondo insulator. We also demonstrate that a correlated 
metallic state with characteristic properties emerges in 
the paramagnetic phase slightly above the AF transition 
temperature, when the Zeeman energy is nearly equal to 
the spin gap. 

This paper is organized as follows. In the next section, 
we briefly mention the model and method, and show the 
phase diagram obtained by DMFT. In Sec. Ill, we explain 
the detail of the field-induced AF transitions, and discuss 
some remarkable properties in the paramagnetic phase. 
Brief summary is given in Sec. IV. 



II. MODEL AND PHASE DIAGRAM 

Let us begin with the periodic Anderson model with 
the Zeeman term, 
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where Ci CT (fi a ) annihilates a conduction (/) electron 
on the ith site with spin a, nf c = a\ a a,ia and Sf z = 
h( n it ~ n i{) w ith a = c, /. Here, t is the nearest- neighbor 
hopping matrix in the conduction band (referred to as the 
c band in the following) . We consider the hyper-cubic lat- 
tice with a bipartite property in infinite dimensions. The 
bare density of states for c electrons is Gaussian with the 
width £*: p (e) = exp[-(e/t*) 2 ]/Vnt* 2 . The magnetic 
field B applied along the z direction is coupled to both 
of the c and / electrons. We assume that the g factors 
of the c and / electrons are the same, and set g^B = 1- 
These simplifications would not qualitatively change our 
conclusions. 

DMFT is a powerful framework to study 
strongly correlated electron systems, e.g., the 
single-band Hubbard model t 24 i 25 ' 26 i 2 7i 28 i 29 i 3n the 
two-band Hubbard model, 31 ! 32 . 33 ! 34 . 35 ! 36 . 37 ! 38 ! 39 the 
PAM, 11 ! 12 ! 14 ! 15 ! 18 ! 20 ! 40 ! 41 ! 42 ! 43 ! 44 ! 45 etc. This method is 
justified in the limit of large spatial dimensions and gives 
a rather good approximation even in three dimensions. 
In DMFT, the proper one-particle self-energy, which 
is independent of the momentum, is obtained via an 
effective impurity model embedded in a self-consistently 
determined medium. The DMFT treatment of the PAM 
is summarized in Appendix A. We solve the effective 
impurity model using the QMC method with the Trotter 
time slices, At = (3/L < 0.25 (/3 = 1/T). The number of 
QMC sweeps is 3.2 x 10 6 in each DMFT iteration loop. 
We use the typical parameters U/t* = 2.0 and V/t* = 0.6 
in the following calculation. For these parameters, the 
ground state at zero field is a paramagnetic Kondo 
insulator, as shown by Jarrell et. alw^ It is convenient 
to express the energy normalized by the zero field value 
of quasi-particle gap, Ao ~ 0.145t*, which is estimated 
from the /-electron spectral function at T/t* = 1/30 
(see Fig. . 

We first give an overview of the phase diagram ob- 
tained in this paper, which is shown in Fig^ At weak 
fields, the system is in the Kondo insulating phase char- 
acterized by the gap formation of quasi-particle states, 
which is smoothly connected to the high-temperature 
paramagnetic phase. Beyond a certain critical field, there 
is a second-order phase transition from the paramagnetic 
phase to the low-temperature AF phase with the decrease 
of temperature. The AF phase possesses a long-range or- 
der in the xy plane, namely, the transverse AF order. 

It is seen that the transition temperature T c takes a 
maximum, when the Zeeman energy is nearly equal to 
the spin gap B ~ 2Aq. Around this field, we find remark- 
able properties in the paramagnetic phase just above T c : 
they are characterized by two crossover fields, Bq and 
Bk- The quasi-particle gap is closed for fields higher 
than Bq. At Bk , the magnetization of the c electrons 
changes its sign. The region surrounded by T c , Bq and 
Bk may be referred to as a 'Kondo metal', where the 
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FIG. 1: Phase diagram of the PAM on the hypercubic lattice 
in infinite dimensions. The parameters are chosen as U = 2.0 
and V — 0.6. The second-order AF transition temperature 
Tc is indicated by open circles. Two fields Bg and Bk char- 
acterize the crossover behavior related to the gap closing and 
the change in the magnetization process, respectively. The re- 
gion denoted as a Kondo metal possesses interesting physical 
properties different from ordinary heavy electrons, see text. 



quasi-particle gap is closed but the Kondo correlations 
are still dominant to produce the strong renormalization 
effects accompanied by the large mass enhancement. As a 
result, this Kondo-metal region exhibits interesting prop- 
erties coming from the competition of the Kondo corre- 
lations, the AF instability, etc. 



III. 



DETAIL OF THE RESULTS 



We now explain how the above phase diagram is ob- 
tained in the DMFT framework. We first investigate the 
field-induced AF transitions, and then discuss some char- 
acteristic properties in the paramagnetic phase. 



A. Phase transitions 

An AF transition is signaled by the divergence of the 
staggered spin susceptibility Xxx(Q) with Q = [tt, 7r, . . .], 
where the suffix x denotes the direction perpendicular to 
the field. In DMFT, we can calculate this quantity via 
the impurity model, since the irreducible vertex function 
as well as the self-energy is independent of the momen- 
tum in infinite dimensions. 11 The detail of the calculation 
is summarized in Appendix B. 

In Fig. |21 the staggered spin susceptibility Tx*"^(Q) 
is plotted as a function of temperature for different mag- 
netic fields. In the DMFT framework, the susceptibility 
exhibits mean-field behavior, \xx 

(Q) cx [T-TJ- 1 . Thus, 
we evaluate the transition temperature rather accurately 
by extrapolating the inverse susceptibility to zero as a 
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FIG. 2: Temperature dependence of the total staggered sus- 
ceptibility Xxx(Q) with Q — [n,n, ...]. Inset shows the in- 
verse of the susceptibility for / electrons xlx(Q)- Both of 
XxxiQ) an d X^(Q) diverge at the same critical temperature. 
XxxiQ) is wel l fitted with the form l/xL(Q) T — T c , which 
enables us to determine the transition temperature T c by a 
proper extrapolation procedure. 

function of T, as shown in the inset of Fig. The tran- 
sition temperature T c thus obtained for various values of 
B is plotted in the phase diagram shown in Fig. ^ The 
magnetic field triggers a phase transition from the para- 
magnetic Kondo insulator to the transverse AF ordered 
state. The critical field between the Kondo insulator and 
the AF insulator is estimated as B ~ 0.56Ao at T = 0. 
In our calculation, it is somewhat difficult to precisely de- 
termine the higher critical field separating the AF phase 
and the paramagnetic metallic phase at T = 0, since the 
signal of antiferromagnetism appears at very low temper- 
ature with increasing B/Aq. Milat et. al. concluded that 
the AF state should persist up to the fully polarized state 
because of perfect nesting. 22 Our numerical results sup- 
port their conclusion, although such a transition would 
occur at extremely low temperatures in the high field 
regime. It is instructive to note that the transition tem- 
perature T c takes a maximum value T c ~ 0.22Aq around 
the field B ~ 1.93Ao, which is almost the same as the 
spectral gap, 2Ao. We will see that this value of field is 
close to two typical crossover fields discussed below. 



B. Heavy- fermion behavior in paramagnetic phase 

We obtain dynamical quantities by applying the max- 
imum entropy method (MEM) to the imaginary-time 
QMC data computed with DMFT. In Fig. we show the 
one-particle spectral functions for the c and / electrons, 
A c ^i(uj) and Af^i(uj), at T = 0.23A for different mag- 
netic fields. Note that the system is in the paramagnetic 
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FIG. 3: Spectral functions of / (left) and c (right) electrons 
at T = 0.23Ao for different values of B/Ao- The energy scale 
Ao is determined by the distance between zero frequency and 
the frequency which takes the half value of the first peak in 
the /-electron spectral function at zero field. 



phase for any magnetic fields at this temperature. One 
can see in this figure how the magnetic filed closes the 
quasi-particle gap and then induces correlated metallic 
states. At B = 0, the well-defined insulating gap ex- 
ists both for / and c electrons. As the field is increased 
(B/Ao = 0.83), the peak structure existing beside the 
gap is slightly shifted in the presence of the Zeeman split- 
ting. For B/Ao ~ 2, the spectral gap almost disappears. 
A remarkable point in this case is that / electrons get an- 
other sharp peak structure around u> = 0, as seen for the 
third panel in Fig. [3J This implies that strong Kondo 
correlations still persist even in such fields, driving the 
system to a correlated metallic state (Kondo metal in 
Fig. ^| immediately after the quasi-particle gap is closed 
by the field. If the field is further increased, such a sharp 
quasi-particle peak disappears, and the system gradually 
changes to an ordinary paramagnetic metal with the spin 
polarization. 

In order to see the above characteristics clearly, we 
show the /-electron spectral weight at the Fermi level, 
A f (0) = [A fl (0) + A h (0)]/2, in Fig. H (a) as a function 
of the magnetic field at different temperatures. It may 
be legitimate to define the crossover field Bq character- 
izing the gap-closing by the field where ^4/(0) takes its 
maximum as a function of B. The crossover field thus es- 
timated, Bq ~ 2Aq, is consistent with the known results 
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FIG. 4: (a) The /-electron spectral function at the Fermi level 
Af(0) as a function of B/Ao, which is used to determine Be; 
(b) the derivative of the /-electron self-energy Z(B)^ 1 at T = 
0.23Ao normalized by its zero- field value Z(0) -1 (~ 2). Here 
Z is computed numerically by the formula Z = ImE(wo)/wo 
with ljq — tyT. 



that the spin gap is nearly equal to the quasi-particle gap 
in infinite dimensions piL^ii It is seen that Bg is almost 
independent of the temperature in the range shown in 
Figffl As a reference, in Fig. 0] (b) we show the nu- 
merical derivative of the /-electron self-energy, Z, which 
would correspond to the inverse of the mass-enhancement 
factor at low temperatures. It is seen that Z~ l indeed 
takes its maximum around the field Bq, which seems con- 
sistent with the development of the sharp quasi-particle 
peak in Fig. [3J However, one notices that the mass en- 
hancement (i.e. development of the sharp peak) observed 
around B/Ao ~ 2 in Fig. [31 is larger than that in Fig. 
0Jb). This difference comes from our numerical differ- 
entiation procedure with discrete frequency in Fig. HJb), 
where we have naively used the formula in the Matsubara 
form, Z = ImS(LUo) / loq with ujq = ttT. This formula may 
not pick up the sharp structure of the spectrum around 
B/Ao ~ 2 in Fig. [3J Therefore, the actual mass enhance- 
ment should be somewhat larger than that displayed in 
Fig. Eft). 

We further calculate the dynamical spin-correlation 
function to directly observe the magnetic-field effects on 
spin excitations. The imaginary-time spin-correlation 
function X ™(r), defined as X M = (T T 5™ (t) (0)) 
with m = c, / and a = z, x, is calculated by the QMC 
method and is analytically continued to the real fre- 
quency with MEM. In Fig. we show the field depen- 
dence of Imx™(u>) at a given temperature, T/Aq = 0.23, 
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FIG. 5: Dynamical spin-correlation function, Imx™ (w), for 
/ (left panel) and c (right panel) electrons at T/Ao = 0.23 
for different values of B/Aq. Solid (broken) lines represent 
the z- (x-) component of the spin correlation function, where 
the z axis is parallel to the applied field. 



slightly larger than the maximum AF transition temper- 
ature. 

At zero field, there are two characteristic features in 
the dynamical spin-correlation functions. One is a low- 
energy peak (shoulder) structure in Imx^(w) (Imx„(w)) 
due to spin-triplet excitations, which determine the size 
of the spin gap. The other is a broad continuum in 
the high-energy region due to the dispersion of c elec- 
trons, which features a higher-energy hump in Imx^(w). 
When a magnetic field is introduced, the spin-triplet ex- 
citations split into three distinct excitations specified by 
Sl ot . At B/A = 0.83, the magnetic field little affects 
the shape of Imxl' c (<jj), whereas it splits the low-energy 
peak of lmxl' c (uj) into two sub-peaks. The results are 
naturally understood, because the triplet excitations ob- 
served in lmxl' c (oj) (Imx£ c (w)) are those with Sl ot — 
(Sl ot — ±1): the latter causes the splitting in the spin 
excitation spectrum. If the field is further increased, 
the spin gap disappears around B / Aq ~ 2, making the 
Kondo singlet state unstable. Note that around this 
field, Imxl(uj) develops another hump structure around 
u) = 0, which is consistent with the fact that heavy- 
fermion states with the enhanced mass appear in the one- 
particle spectrum (Fig. Another noticeable point is 
that ImXzi^) does not increase its weight around lu = 0, 
although Im^(w) as well as lmx( x (uj) already has a 
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FIG. 6: The magnetization process of / (a) and c (b) elec- 
trons. Mf and M c are defined by Mj = (n/j — n/j.) and 
M c = (n cT - n c i ) . 



reasonable weight (see data at B/Aq = 2.07). Further- 
more, lmx c z (w) eventually develops a shoulder-like struc- 
ture around lu = at B/Aq = 3.31, where the heavy- 
fcrmion behavior in lmxl(oj) is already suppressed. In 
this way, the z-component spin excitations of c electrons 
exhibit unusual behavior around the field B/Aq — 2^3. 

To make the above point clear, we turn to another in- 
teresting crossover behavior denoted by the line Bk (see 
FigU in the paramagnetic phase above T c . This line is 
determined by the change in the magnetization process 
for c electrons. Shown in Figs. El (a) and (b) are the 
magnetization of / and c electrons at several different 
temperatures. It is found that the magnetization of / 
electrons is always positive, while that of c electrons is 
negative for weak fields and gets positive beyond a cer- 
tain magnetic field. We define Bk by the magnetic field 
at which the magnetization of c electrons changes its sign. 
This field Bk roughly corresponds to the strength of the 
local singlet formation. At weak fields, the Kondo effect 
is still dominant to align the spins of c and / electrons in 
the opposite directions, which may give rise to the neg- 
ative magnetization in c electrons. Such nonmonotonic 
behavior is also seen in the spectral function of c elec- 
trons in Fig. |3| The spectral function of c electrons with 
up spin is shifted to the low-energy side, while its weight 
is considerably suppressed. These results are also consis- 
tent with what we have observed in the small-w behavior 
of ImXcH for B/A = 2 - 3 in Fig. 

Summarizing the results found in the normal phase, we 
can say that around the field where T c takes a maximum, 



various correlation effects such as the Kondo-gap forma- 
tion, the AF instability, etc. compete with each other, 
yielding characteristic properties in correlated electrons. 
Such features are different from ordinary heavy fermions 
without magnetic fields. 



IV. SUMMARY 

We have investigated the effects of the magnetic field 
on the Kondo insulator by applying DMFT to the PAM 
at finite temperatures. There are several characteristic 
features in the phase diagram. At low temperatures, we 
have observed the quantum phase transition from the 
Kondo insulator to the transverse AF phase, in accor- 
dance with the results obtained in the two-dimensional 
case. We have also found some remarkable properties in 
the paramagnetic regime. Namely, in the paramagnetic 
phase close to the maximum AF transition temperature, 
heavy-fermion states with the enhanced mass are real- 
ized, where the competition of the magnetic field and 
the Kondo correlations results in characteristic features 
in the one-particle spectrum, the dynamical spin suscep- 
tibility and the magnetization process. It may be inter- 
esting to experimentally investigate static as well as dy- 
namical properties in this region. Such investigation may 
provide a unique example of field-induced heavy-fermion 
states without phase transitions. 
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APPENDIX A: FORMULATION OF DMFT 

In this Appendix, we summarize the DMFT formula- 
tion of the PAMi 11 ! 12 ! 13 In DMFT, the original lattice 
model is mapped onto an effective impunity model with 
the effective medium determined self-consistently. This 
medium is described in terms of a cavity Green's func- 
tion Q(iuj n ), which represents the motion of correlated 
electrons on all sites except an impurity site. The local 
Green's function G(iui n ) is then expressed in terms of 
the momentum-independent self-energy E(ia>„) and the 
cavity Green's function Q(iuj n ) as 

G°{iuj n ) = [Q a {iu n )- X - VilLUn)]- 1 , (Al) 
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where the Green's functions and the self-energy are ex- 
pressed in the matrix form, 



G%(iuj n ) G° }f {iuj n ) 



E?(iw„) 



(A2) 
(A3) 
(A4) 



On the other hand, the local Green's function G(ioj n ) 
can be written as 



G°(iu n ) 



de 



go(g) 



£?(*w„) V 



with 



£° {iu n , e) = iu n + -<Tgp B B - e, 



(A5) 



(A6) 



C/(iw n ) = ioj n - — + ^ag/j, B B - S/ (iw n ) ,(A7) 



(A8) 



By solving equations IjAlj) and l|A5|l self-consistently, we 
can discuss the effects of electron correlations on local 
properties in the PAM. In this paper, we obtain the 
imaginary-time Green's function G" T (r) by solving the 
impurity problem with the QMC method^ 



APPENDIX B: MAGNETIC SUSCEPTIBILITY 

To discuss the instability to an AF ordered phase, it 
is necessary to obtain the momentum-dependent suscep- 
tibility Xxx(<l) in the PAM. Here, we explain how this 
quantity is formulated in the framework of DMFT>ii*H 

We first consider the local two-particle propagator in 
the effective impurity system, 



X+-i (n,T 2 ,T 3 ,r 4 ) 



T T f^( Tl )f U (r 2 )4(r 3 )/. iT (r 4 )),(Bl) 



X- +t (ti,t 2 ,T3,t 4 ) 



Tr/i(ri)/i T {n)flr(n)fn (r 4 )).(B2) 



We obtain the vertex function T(iui = 0; iu) mi iu) n ) by 
solving the Bethe-Salpeter equation, 



% A% 1 A% At i 



(B3) 



where x{° 1S the non-interacting part of the local two- 
particle propagator, 

x{° (m = 0; iuj m , iw n ) = — ( iuj n) G l ff (iio n ) . 



(B4) 



By using the vertex function T, we obtain the two-particle 
propagator for correlated /-electrons in the lattice sys- 
tem, 



yS = y /0 , y/0 ry / 
Aq Aq 1 Aq Aqi 



(B5) 



where x q ° is the bare two-particle propagator in the lat- 
tice system given by, 



X+°- q = 0; iu m , ioj n ) 



N/3 



^G%(iu n )G{^(iu> n ), (B6) 



5. 



Nj3 



^G{{(^„)G^ qT (^ n ). (B7) 



Equations l|B6|) and (|B7|) can be rewritten as 
dx 1 dx 2 Aq{xi,X2) 



X-+q(ivi = 0;iuj m ,iuj n ) 



dxidx 2 Aq(xi,x 2 ) 



m i %2 ) 

(B8) 
(B9) 



with the two-particle density of states A q (xi, x 2 ) defined 



as 



Aq(xi,X 2 ) = -^^2${ x l ~ £k)$( x 2 - £k+q)- (BIO) 



If we focus on the center or the corner of the Brillouin 
zone, this function is reduced to the simpler form, 



Aq(xi,a;2) = 



_ J 6(xi - x 2 )po(xi), q = [0, 0, . 



S(x 1 +x 2 )po(x 1 ), q = [it, tt, . . .] 



(Bll) 



Using the above equations, we can calculate the 
momentum-dependent two-particle propagator Xqfal = 
0; iuj m , iuj n ). We then obtain the static susceptibility, 

x/(q) = p^Xqim = 0;iw m ,iw n ). (B12) 



The total two-particle propagator in the lattice system, 
defined as x q ot = X q + Xq + X q c + X q > can be expressed 
by that for correlated /-electrons xL as 



y tot = c0_ c/0/ /OA- 1 /c0 
Aq Aq Aq \Aq ) Aq 



i + x q /0 (x q T 



Xq 



i+(x q °) xr 



(B13) 



7 



where I the unit matrix, and X q °(^z = 0;iu> m ,iuj n ) ancl 
Xq^°(*^i = 0;iw m ,iw„) are the c-c and c-f elements of 
the bare two-particle propagator, which are obtained in 



the same way as x q °(^z — 0;icj m ,iu) n ), We then obtain 
the total static-susceptibility x tot (q), summing over the 
Matsubara frequency in Xq°*(^( — 0; ** w m> i^n)- 
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